function [bgRegionStd,bgRegionMean]=ODimageMeanStd(odimage)
    %cal standard devation of background
    bgROI=[667 191 300 600];
    bgRegion=odimage(bgROI(2):bgROI(2)+bgROI(4),bgROI(1):bgROI(1)+bgROI(3));
%     odimage(bgROI(2):bgROI(2)+bgROI(4),bgROI(1):bgROI(1)+bgROI(3))=10;
%     figure(10);imagesc(odimage);
    bgRegionStd=std(std(bgRegion));
    bgRegionMean=mean(mean(bgRegion));
    rectangle('position',[bgROI(1) bgROI(2) bgROI(3) bgROI(4)],'Curvature',[1,0],  'EdgeColor',[1 0 0]);